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A Bose-Einstein condensate in an optical lattice exhibits parametric resonances when the intensity 
of the lattice is periodically modulated in time. These resonances correspond to an exponential 
growth of the population of counter-propagating Bogoliubov excitations. A suitable linearization of 
the Gross-Pitaevskii (GP) equation is used to calculate the stability diagram and the growth rates of 
the unstable modes. The results agree with the ones extracted from time-dependent GP simulations, 
supporting our previous claim (M. Kramer et al, Phys. Rev. A (2005) in press) concerning the 
key role of parametric resonances in the response observed by Stoferle et al. (Phys. Rev. Lett. 
92, 130403 (2004)) in the superfluid regime. The role of the seed excitations required to trigger the 
parametric amplification is discussed. The possible amplification of the quantum fluctuations present 
in the quasiparticle vacuum, beyond GP theory, is also addressed, finding interesting analogies with 
similar processes in nonlinear quantum optics and with the dynamic Casimir effect. Our results can 
be used in exploiting parametric instabilities for the purpose of spectroscopy, selective amplification 
of a particular excitation mode and for establishing a new type of thermometry. 

PACS numbers: 03.75.Kk, 03.75.Lm 



I. INTRODUCTION 

A parametric resonance corresponds to the exponen- 
tial growth of certain modes of a system induced by 
the periodic variation of a parameter 0, 0- It is 
a very general phenomenon occurring in classical os- 
cillators, in nonlinear optics, in systems governed by 
Non-Linear Schrodinger equations, and in Hamiltonian 
chaotic systems. Parametric resonances can also occur 
in Bose-Einstein condensates made of ultracold atomic 
gases. From the theoretical side, they have already 
been investigated in 3D 0. [H and 2D 0, Q harmon- 
ically trapped condensates, for condensates in oscillat- 
ing double-well potentials [8|, la, |lfj and in deep optical 
lattices These investigations are based on suitable 
effective Hamiltonians in which the interaction between 
the atoms enters through a mean-field term proportional 
to the condensate density. Most of them make use of 
the Gross-Pitaevskii (GP) equation, which has the form 
of a Non-Linear Schrodinger equation and gives a very 
accurate description of dilute condensates [12J, Il3j | . 

In Ref. 01 we have recently shown that parametric 
resonances can be the origin of the large and broad re- 
sponse observed by Esslinger and co-workers 0, 0, 
in the superfluid phase of a condensate in an optical lat- 
tice. Our claim was based on the results of numerical sim- 
ulations, namely the integration of the time-dependent 
GP equation. In the present work we perform a more 
systematic analysis, by investigating the behavior of Bo- 
goliubov quasiparticles subject to a periodic modulation 
of the lattice. The purpose is to gain a deeper insight 
into the process of parametric resonances, in particular 



into the conditions for the occurrence of the instability, 
the growth rates associated with the parametric ampli- 
fication of an excitation and the role of the initial fluc- 
tuations (seed excitations) which trigger the onset of the 
amplification. 

The basic theory is presented in section|n] We consider 
an infinite condensate at zero temperature, subject to a 
transverse harmonic potential and to a ID optical lat- 
tice along the z-direction, with TV atoms per lattice site. 
The effect of the optical lattice can be expressed through 
the periodic potential V(z,t) = s (t) E r sin 2 (qsz), where 
qB = ir/d is the Bragg wave vector, d is the lattice spac- 
ing and s is the lattice depth in units of the recoil en- 
ergy Er — Tiujr — h 2 q B /2m. The transverse harmonic 
confinement, of frequency uj±, is taken to be sufficiently 
strong to inhibit excitations involving the radial degrees 
of freedom on the energy scale we are concerned with. 
Under this assumption, the order parameter can be fac- 
torized into a radial Gaussian, having constant width 
aj_ = [ft,/ '(mcjjj] 1 / 2 , times a z and i-dependent order 
parameter, ^>(z,t). The latter obeys an effective ID 
GP equation which includes the atom-atom interaction 
through a coupling constant g. We first consider the case 
of a static optical lattice (s(t) — so) and we find the Bo- 
goliubov spectrum as a result of the linearization of the 
GP equation with regard to a small change, 5' l S(z,t), of 
the condensate wave function. Then we study the evolu- 
tion of the system under a modulation of lattice depth in 
the form s(t) — sq[1 + Asm(ilt)]. A stability analysis is 
done in order to find the parametrically unstable modes 
and calculate their growth rate for given values of the 
parameters so, A and gn, where n = N/d is the average 
linear density. We use a Bogoliubov quasiparticle projec- 
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tion method |21j, |22j in order to provide a deeper insight 
into the mechanisms of instability. 

When A is not too large, the main mechanism which 
drives a parametric resonance is a coupling between pairs 
of counter-propagating Bogoliubov axial excitations of 
frequencies u>jk and u>j'k, where k is the wave vector and j 
and f are Bloch band indexes. The multi-mode stability 
analysis then reduces to a two-mode approximation. The 
resonance condition becomes — ujjk+ujj'k, which yields 
O = 2cjofc for axial phonons in the lowest Bloch band. 
This approximation is discussed in section ITTT1 where we 
show that the parameters so, A and gn can be combined 
in a single parameter 7 = (l/2)soAD,T^_, which con- 
tains all relevant information about the stability diagram 

and the growth rate. The quantity r+_ depends on so 
and k and can be calculated by solving the Bogoliubov 
equations for different lattice depths. It describes the 
coupling between Bogoliubov modes belonging to the j 
and j' bands and having opposite quasi-momenta. This 
coupling is mediated by the time-dependence of the lat- 
tice potential. In the tight binding regime, i.e. when the 
lattice depth is large enough to split up the condensate 
in almost independent condensates at each lattice site, a 
semi-analytic result can be derived for the quantity 7 de- 
scribing the parametric resonance of modes of the lowest 
Bogoliubov band. In this case, in fact, the quantity r+°_ 
can be simply expressed through the s-dependence of the 
effective mass and the compressibility of the condensate. 
This regime is discussed in section llVl 

The numerical results of the stability analysis are 
shown in section together with those of the two- mode 
and tight- binding approximations. The two- mode ap- 
proximation turns out to be very accurate in the whole 
range of parameters here considered. At large So the 
tight binding expressions work well. The approach to 
this regime is faster when the mean-field interaction is 
small (i.e., small values of gn). 

In section IVTl the results obtained in the previous sec- 
tions by linearizing the GP equation for an infinite con- 
densate are compared with those obtained from the nu- 
merical integration of the time dependent GP equation. 
We perform GP simulations for an infinite condensate, 
as well as a trapped condensate similar to the ones of 
the experiments of Refs. [l5llT^.lr7| and of our previous 
calculations of Ref. [l4|. The agreement is good. This 
confirms the overall picture for the key role played by 
the parametric instability of Bogoliubov excitations in 
the observed large and broad response of the condensate 
to the modulation of the lattice depth. These results also 
support the possibility to use parametric resonances as a 
tool for spectroscopic studies of Bogoliubov excitations, 
for a selective amplification of certain modes and for the 
characterization of thermal and/or quantum fluctuations 
in trapped condensates. A brief discussion on how para- 
metric resonances can give information about thermal 
and quantum fluctuations is given in section lYm 



II. BOGOLIUBOV EXCITATIONS IN AN 
AMPLITUDE MODULATED LATTICE 

A. Bogoliubov spectrum in a static lattice 

Let us consider a dilute condensate made of atoms with 
mass m and s-wave scattering length a, whose order pa- 
rameter obeys the time-dependent GP equation |12l [jjj 



-— V 2 + V + g 3D N\% 2 
Zm 



(1) 



where ^ is normalized to 1, TV is the number of particles 
and gzD = iirh 2 a/m. The potential V is the sum of the 
harmonic trap and the optical lattice. The ID optical 
lattice is taken to be oriented along the z-direction with 
lattice depth s, spacing d, Bragg wave vector qs = 7r /d 
and recoil energy Er = huipi — h 2 q B /2m. The transverse 
harmonic confinement is assumed to be sufficiently strong 
to freeze the radial degrees of freedom, so that the order 
parameter is simply a Gaussian in the transverse direc- 
tion, having constant width aj_ = [H/ '(mwjj] 1 / 2 . Under 
this condition, the 3D GP equation Q reduces to an ef- 
fective ID GP equation for a purely axial order parameter 



2m 



d 2 + sE R sm 2 (q B z) + gnd \^\ 2 



*■ (2) 



Here the normalization condition is J_^i 2 dz\^>\ 2 = 1, 
while nd — N is the number of atoms per lattice site 
and n is the average linear density. The quantity g = 
g 3 £ ) /2Tra 2 L is an effective ID coupling constant accounting 
for the atom-atom interaction. With this choice, Eq. 
also coincides with that of Ref. 01 f° r a uniform con- 
densate in a ID lattice. 

When the condensate is weakly perturbed, the order 
parameter can be written as 



9(z,t) 



— i^t/U 



[* o (a0 + **(*,*)] 



(3) 



where ^0 is the groundstate solution of the stationary 
GP equation 



2m 



d 2 + sE R sin 2 {q B z) + gnd |* 



*o = A**o (4) 



and n is the corresponding chemical potential. The small 
variation 5^ can be expressed in terms of Bogoliubov 
excitations, 



5^ 



jk 



■iujjkt 



c jk v jk 



(z)e 



iujjkt 



(5) 



where Cjk are complex coefficients. Due to the periodicity 
of the external potential, the Bogoliubov quasiparticle 
amplitudes Ujk and Vjk have the form of Bloch functions 
Ujk = exp(ikz)iijk and Vjk — exp(ikz)vjk, where Ujk 
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and Vjk are periodic with period d. The excitations are 
labelled by their quasimomentum hk and the Bloch band 
index j. 

By inserting © into 10 and linearizing with respect 
to <5\& one finds that the amplitudes Ujk and Vjk are so- 
lutions of the eigenvalue problem 

where the matrix £(s, fc) is given by 

M*,*)- {-gncMf -H(s,k) ) (7) 
and we have defined 

h 2 

H(s,k) = - — {d z +ik) 2 + sE R sm 2 (q B z) 
Am 

+ 2gnd |^o| 2 — /it . (8) 

Moreover, for any given k the solutions Ujk and Sjfc obey 
the ortho-normalization conditions 

fd/2 

/ dz(u* jk Uj, k -v* jk v rk ) = Sjj> , (9) 

J-d/2 
fd/2 

/ dz^jkVj^-k-uj'-kVjk) = 0. (10) 

J-d/2 

Equations © can be solved numerically to get the am- 
plitudes Ujk, Vjk and the Bogoliubov band spectrum tujjk 
for given values of the lattice depth s — s and interaction 
parameter gn. In Fig. ^ we plot the lowest three bands as 
obtained for sq = 4 and gn = 0.72Er. The latter value is 
chosen as the axial average of the interaction parameter 
in the experimental setting of |Tfij |. The lowest band ex- 
hibits a linear dispersion for k — > as expected for long 
wavelength Bogoliubov phonons traveling along z at the 
speed of sound. This type of spectra have been already 
discussed in the recent literature 0, 0, |2(| . Here we 
use the results of Fig. ^ as a basis for the subsequent 
discussion of parametric resonances. 



B. Modulation of the lattice depth 

Let us consider a periodic modulation of the lattice 
depth in the form 

s(t) = s [1 + A sin(fii)] , (11) 

with A -C 1. The dynamics induced by this modulation 
can still be described by Eqs. @ and ©. To this purpose 
we assume that \&o is, at any time t, the stationary solu- 
tion of (0} for s = s(t). Both ^ and ^ are now functions 
of t via their dependence on s(t). Then we insert again 
(J3J into J2J) take s as in (|llf> . and look for an equation 
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FIG. 1: Lowest bands in the Bogoliubov spectrum of a cylin- 
drical condensate in a periodic lattice, obeying Eq. with 
s = so = 4 and gn — 0.72Er. The range of k corresponds to 
half of the first Brillouin zone. We use the shorthand nota- 
tion u>o, wi and u>2 to indicate the Bogoliubov bands whose 
dispersion relation is u>ok, Wit and oj^k- 



of motion for the small part Sty . Up to linear terms in 
(5 'J, one gets 

idttftf = [H(s,k) + 2gnd |* | 2 - 

+ gnd^ 2 8^* + s[id s + (d,n)t/h]9 ■ (12) 

The inhomogeneous term, proportional to s, behaves like 
a source of excitations in the linear response regime. It 
creates excitations out of the condensate when fi is res- 
onant with a Bogoliubov frequency ujjk at k = 0. For 
symmetry reasons, the first resonance occurs at the bot- 
tom of the 0J2 band. In the case of Fig. ^ this amounts 
to fl = uj2(k — 0) ~ bAEn/h. Similar resonances can be 
found at larger frequencies. If A is small and the modula- 
tion time is long, this type of excitation processes only oc- 
curs in narrow and well separated intervals of f2. Except 
in these intervals, the inhomogeneous term in Eq. I|12|) 
can be safely ignored. In the range of of the two low- 
est Bogoliubov bands loq and to>i, for instance, this ap- 
proximation is expected to work well, as confirmed also 
by the comparison with GP simulations (see section IVT1 
below). In the same range, also the assumption that 
adiabatically follows the lattice modulation is expected 
to be valid. Finally, the periodicity of the system allows 
one to write the unknown function S^(z,t) in the form 
of a Bloch wave expansion 

<y*(z,t) = Y,a+k{z,ty kz + ~at k { Zl t)e~ lkz , (13) 

fe>0 

where the coefficients a k (z,t) are periodic in z with pe- 
riod d. Neglecting the last term in (|12|) and inserting the 
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ansatz <|13|) . one finally gets 



C. Stability analysis 



ihdt 



a+k 
a-k 



= C(s,k) 



a+k 
a-k 



(14) 



This is the equation of motion for 5ty. It contains the 
external parameter s(t), which varies periodically in time, 
and hence it can exhibit parametric instabilities. We 
perform the stability analysis of this equation using the 
method explained in section III CI below. 

In order to give a more transparent view of the pro- 
cesses involved in Eg . II 1411, one can use a suitable basis 
of Bogoliubov states |2ll |22| . Let us write 



a+k 
a~k 



E 

3 



Ujk 

v jk 



(15) 



where iijk, Vjk denote the solutions of (jSJ at the lat- 
tice depth s(t). This corresponds to the projection of 
Sty, at time i, onto the instantaneous Bogoliubov modes 
at lattice depth s(t). Inserting (|15fl into l|14l) and using 
(|§|)- (I1L)|> . one finds the following equations for the time- 
dependent coefficients Cjk- 



id t Cjk = iOj k c ]k -is) j \ T J l + c rk + T+- 



(16) 



id t c* 



where 



j' 



Cfk + 



V 33 ( 
1 ++ c 



-p33 
1 +-f 



d/2 

dz (u* 



jk 



■ d s u 



1 +- 



-d/2 
d/2 

dz (u* k d s v*, 

-d/2 



j'k 



v* k d s u*. 



-k 



(17) 



(18) 



(19) 



In deriving (|16l) - (|17fl we have used the relations 



and Vi 



ijk - 

which follows from the fact that 



~j,-k allu V jk — V j,-k' 

the groundstate wavefunction can be chosen real (zero 
current). 

Eqs. Hltjfl - I|17|) reveal that the time-dependence of the 
lattice depth (s ^ 0) brings about a coupling between 
Bogoliubov excitations with the same In particular, 

the quantity T 3 P, accounts for the scattering of excita- 
tions from the perturbation, that is a process in which 
an excitation in the band j is scattered into the band j' 
without changing its quasimomentum k. This happens 
resonantly when f2 is equal to the frequency difference 
between two bands at the same k. One can easily prove 
that this process does not vary the total number of exci- 
tations [23J. Conversely, the quantity r+_ describes the 
coupling between counter-propagating excitations. As we 
will show in the following, this quantity causes an expo- 
nential growth of the number of excitations and is thus 
responsible for the parametric instability. It is worth 
noticing that T 3 _^_ = if gn = and hence a noninter- 
acting gas, obeying the Schrodinger equation, does not 
exhibit parametric resonances. 



In order to find the parametrically unstable solutions of 
the linear equation (|14|) we follow the procedure discussed 
in Ref. 24] . The dynamics is determined by the operator 
£, which changes in time via its dependence on s(t). The 
lattice modulation is periodic with period T — 2-k/Q. 
Therefore the operator C has the same periodicity and 
can be expressed by means of the Fourier representation 



+ OC 



£(t)= £ Mg 



(20) 



v— — oo 



where v is integer and the operators do not depend 
on time. Thanks to Floquet theorem [25|, any solution 
of (|14|) can be written in the form 



a+k(t) 
d-k(t) 



a'+ k {t) 
a'- k (t) ) ' 



(21) 



where A is a complex constant and the functions a'±k{t) 
are periodic of period T . This periodicity allows one to 
expand a'±k in the Fourier series 



+oo 



a'±k(t)= 



(22) 



v— — oo 



Finally, by using Eqs. lf2U|l - l(^|) one can rewrite Eq. l(Ti|) 
in the form 



iX 



~0) \ 

a -k J 




From the solutions of this eigenvalue problem one obtains 
the complex eigenvalues A as a function of k and Q. In 
practice, the solution can be found by truncating the 
series l|22fl at a certain v so that the problem is reduced 
to a matrix diagonalization. In our case, we obtain a very 
good convergence by cutting at v of the order of 10. 

We then define the growth rate 7 of a Bogoliubov mode 
as the real part of A: 



7 = Re(A) 



(24) 



This quantity fixes the character of the corresponding 
eigenfunction. If 7 < then |ci±fc(£)| is always bound 
and we say that the mode ±k is parametrically stable. 
On the contrary, if 7 > the corresponding eigenfunction 
describes an unstable mode whose population exponen- 
tially increases in time. For a given modulation frequency 
we say that the condensate is parametrically stable if 
all eigenfunctions of are stable. By using symmetry 
properties analogous to the ortho-normality relations © 
and IjlOl) , one can prove that the condensate is paramet- 
rically stable if and only if 7 = for all modes. The 
results of this stability analysis will be given in section 
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III. TWO-MODE APPROXIMATION 



with initial conditions 



According to Eqs. I|lfl |l -I|17 [l . the evolution of a certain 
mode is, in principle, coupled to the modes of all bands 
with the same |fc|. However, one might expect that each 
resonance essentially arises from the coupling of a pair of 
counter-propagating modes, the influence of all the others 
being negligible. In our formalism, this is equivalent to 
replacing the sum over j' in Eqs. (|TB|l - l(T7|) with a single 
term. The problem is thus reduced to the solution of the 
coupled equations 



id. 



Cjk 

ft 

k 

- isT J ? _i 




k of the zt/c-modes in the 
is the same 



for the coefficients Cjk and c 
band j and j' respectively. The coupling T J P 
already given in i|19|) . 

The modulation l|ll|) of the lattice explicitly enters 
the two-mode equation (|25[l through the term s — 
Aso^l cos(fii) in the matrix elements. An additional 
implicit dependence arises from the s-dependence of 

In order to investigate the 



IjJjk, Wj'k> 



■ ++ 



and r 



regime of small modulation amplitude, one can linearize 
Eq. (|25|l with respect to A, keeping only the leading order 
of each matrix element. One obtains 



id t 



c jk 

r* 

C j'-k 

ujjk —i2 1 J cos(ttt) 
-i2jcos(Q,t) —ujj>k 



Cjk 

c *r-k 



where 



and 



-^boO(lf_) s=so 



Uljk — {uJjk)s=s 



(26) 



(27) 



(28) 



A further simplification is obtained by neglecting the cou- 
pling between resonances at ±f2, which is justified when 
the modulation time is much longer than f2 — 1 . This 
allows one to replace 2cos(f2i) with exp(ifii), so that 
Eq. yields the two coupled equations 



dtfjk = -7e l H'*+^-^/* 



dtff 



h 



(29) 
(30) 



where we have introduced the new variable fjk = 
e^^Cjk- By inserting (|3U[I in the time derivative of (|29|l . 
one gets 

[d 2 t - i(u jk + Uj, k - (l)8t - |7| 2 ] fjk{t) = (31) 



/ Jfc (0) = c jk (0) ; (d t f jk ) t=Q = -74,_ fc (Q) ■ (32) 

One can easily see that the solutions of Eq. (|31() exhibit 
an exponential growth within the unstable region \ujj k + 
uij'k — f2| < 2|7|, where the growth rate is found to be 

7 =i[4|7| 2 -K fc +^-fi) 2 ] 1/2 . (33) 

This also implies that the maximal growth occurs at the 
rate 

7max - |7l (34) 

when the resonance condition 

= Ujk + LUj'k (35) 

is met. The width of the unstable region is 2| 7 |. By 
recalling definition Ij27(l , one sees that both the maximum 
rate and the width of the unstable region are proportional 
to the modulation amplitude A. We notice that, when 
the coupling occurs between counter-propagating modes 
in the same Bloch band (J = j'), the resonance condition 
is simply £1 = 2u>j k . In this case, by using the properties 
(|9"|)- (llL)ll . it is also possible to show that Im[Ti?_] = and 

If + = o. 

These are the main results of the two-mode approxi- 
mation. We note that the evaluation of the growth rate 
B33JI only requires the solution of the Bogoliubov equa- 
tions © at So yielding the frequencies u>jk,u>j'k an d the 
coupling element 

A further advantage of the two-mode approximation 
is that it allows one to characterize the role of the seed 
excitations. Let us concentrate for simplicity on the case 
of modes satisfying the resonant condition f2 = ujjk+ujj'k. 
The solution of Eq. I|31|) with initial conditions l|32|) has 
the form 



where 



f J k(t)=r ] e^ t +T 1 'e- 



V = 2M0)- S V fc (0)e l 
V' = hc jk (0)+4,-k(°y 



(36) 

(37) 
(38) 



and (f> = phase (7). For long time it ^> 7~J X ), the second 
term in l|36l) vanishes and the population of both reso- 
nant modes have the same exponential growth with rate 
7max- The coefficient r\ is related to the initial popula- 
tion of the Bogoliubov modes (seed excitations) through 
Eq. (|37|) . More precisely, the parametric resonance am- 
plifies a quadrature of Bogoliubov excitations which is 
phase locked with the perturbation (2^, |22j , and rj plays 
the role of an effective seed. It is worth stressing that 
exponential growth, i.e., the parametric instability, only 
occurs if the Bogoliubov modes are already populated 
at the initial time (77 7^ 0) . In section IVIII we will dis- 
cuss how this can be related to thermal and/or quantum 
fluctuations. 
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IV. TIGHT BINDING REGIME 



V. RESULTS 



As the lattice is made deeper one enters the tight bind- 
ing regime where the condensate atoms are strongly lo- 
calized near the minima of the lattice potential and only 
next neighbor Wannier functions have nonvanishing over- 
lap. In this regime, one can derive an analytic expression 
for the lowest Bogoliubov band (see for example [l?H ) 



huj = y/e {e + 2n *) 



(39) 



with e = 2<5 sin(fcc?/2) 2 . The tunnelling parameter 5 is 
related to the effective mass through the relation 5 = 
(2/TT 2 )(m/m*)Efi, while k^ 1 — nd/i/dn is the inverse of 
the compressibility. Furthermore, as discussed in [l8| . 
the Bogoliubov amplitudes of the lowest band become 
proportional to stationary Bloch state solutions ipk of 
the GP equation, u ok = U k ip k , v ok = V k (p k with 



U k = 



V k = 



£ + fUjjQ 



e - hoj 



2y/hLO e 



(40) 
(41) 



One can use these results to calculate the coupling ele- 
ment in Eq.O. The result is 



1 +- 



U k d a V k - V k d s Uk • 
Inserting Eqs. I|4U|) - (|41[I . one obtains 



(42) 



r uu_ = 



uj d s e - epd s u)o 
2euj 

k d s 6 — SdgK -1 
2d (e + 2K- 1 ) 



(43) 



In the tight binding regime the coupling T°^_ is hence 
directly related to the lattice depth dependence of com- 
pressibility and tunnelling (effective mass). Expression 
(|43|l allows for a comparison with the two mode approx- 
imation, as we will see in the next section. 

A nice result is obtained by inserting expression i|43|) 
into Eq. (|27fl. using the definition (|34|l and taking the 
so ^ 1 limit. One finds 



Asnil dc 



In 



hd 



with 



4chd ds 



c hd = (to*k) 1/2 , 



(44) 



(45) 



The latter coincides with the Bogoliubov sound velocity 
of the hydrodynamic approach of Ref. [18|, Therefore, in 
the limit of very large so, where the order parameter \&o 
can be assumed to be time-independent and the fluctuat- 
ing part 5^ can be described in terms of hydrodynamic 
phonons with linear dispersion u>o k = c^dk, the growth 
rate of the parametric resonance turns out to be deter- 
mined by the change in sound velocity induced by the 
lattice modulation. 



In this section, we present the results of the stability 
analysis of section [H] together with the ones of the two- 
mode and tight-binding approximations of sections IIIII 
andllVl 

Let us first discuss the limit of a small modulation am- 
plitude, A — > 0. In this limit the results of the stability 
analysis of Eq. I|14|l are found to exactly coincide with the 
ones of the two-mode approximation. The parametrically 
unstable regions in the (f2, k) plane have a vanishingly 
small width, falling onto the lines where the resonance 
condition (|35|) is fulfilled. The lowest lines are plotted in 
Fig. [5J The different line styles indicate different pairs of 
resonant modes. By direct comparison with the spectrum 
in Fig. ^ one sees that the first resonance is encountered 
at twice the frequency of the lowest Bogoliubov band, 
il = 2cjofc (solid line). The next resonance is found when 
fl is the sum of the frequencies of the lowest and the next 
Bogoliubov band, fl = ujQ k + uj\ k (dashed line). Further 
resonances occur at twice the frequencies of the first ex- 
cited band, = 2ui k (dash-dotted line), at the sum of 
frequencies of the lowest and the second excited bands, 

= uj{) k + uj2 k (dotted line), and so on. One or more 
parametric excitations occur for any value of f2, except 
within the gap between 2a>o and {loq + u)\) resonances at 
zone boundary. 
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FIG. 2: Position of parametric resonances for so = 4, gn = 
0.72Er and for a small modulation amplitude, A — + 0. See 
FigQfor the Bogoliubov spectrum ujj k with the same param- 
eters. 

For small A, both the growth rate of the unstable 
modes and the width of the instability regions are found 
to be proportional to A. In Fig. [3] we plot the growth 
rate on resonance, 7 = 7 ma x> divided by A, for the same 
resonances of Fig. [5] Solid lines are the numerical re- 
sults of the stability analysis of section III Gl (see defini- 
tion J23J) for A < 0.01, while empty circles represent 
the results of the two-mode approximation (see Eqs. I|27|) 
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and JHHO- The very good agreement between the exact 
multi-mode calculations and the two-mode results shows 
that, for small A, the essential mechanism underlying 
the parametric process is fully captured by the two-mode 
approximation. In Fig. [2] the quantity j/A takes values 
ranging from to typically ~ 0.2cjr, which implies that, 
on resonance, parametric amplification occurs on a time 
scale much longer than ui^ 1 . 
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FIG. 4: Stability diagram in the (k, f2)-plane obtained with 
the procedure of section III CI for A — 0.1, so = 4 and 
gn = 0.72Er. The lightest level of grey in the grey-scale 
represents all modes having 10~ 10 < 7 < 10~ 3 . The main 
branch corresponds to the 2oJo-resonance, whose maximal 
growth rate 7 max grows with k as shown in Fig. The other 
weak branches are higher order parametric resonances. 



B J 



FIG. 3: Growth rates of parametric amplification for so = 
4, gn — Q.72Er in the small A limit. The resonances are the 
same as in Fig|5| Lines: results of the stability analysis of 
Eq. 114H . Different values of A in the range A < 0.01 give the 
same curves for j/A. Circles: Eqs. I|27^ and 1)34^ of the two- 
mode approximation, with T^f_ given in 1191 and calculated 
by solving the Bogoliubov equations (JSJ at so = 4. 

Within the two-mode approximation, two excitations 
of frequency u>ik and ujjk are unstable when \u>ik + u>jk — 
Q\ < 2|t|, where |t"| = 7max is the maximal growth rate 
on resonance. The values of the maximal growth rate 
plotted in Fig. [3] thus indicate that, at small values of 
A, one is generally dealing with very narrow resonances. 
For the parameters used in Figs. |21 and |21 (sq = 4 and 
gn = 0.72Er) and for A < 0.01, the width is less than 
1% of the resonance frequency. 

The overall qualitative picture is preserved when the 
modulation amplitude A is tuned to experimentally rel- 
evant values. An example is given in Fig. ^ where we 
plot the stability diagram in the energy range of the 2luq- 
resonance obtained with the stability analysis of section 
III CI for A = 0.1. The main branch remains centered 
around the resonance condition Q — 2ujQk- The shape 
of 7(fc, O) near resonance is very well approximated by 
Eq. 123, which yields 7 = [7^ - {cu ok - ft/2) 2 } 1 / 2 . The 
maximal growth rate 7 max is plotted as a dot-dashed line 
in Fig. El Figure[Hshows that the 2a>o-resonance remains 
narrow even for relatively large A, the finite width be- 
coming noticeable only close to the band edge where it 
reaches about 5% of the resonance frequency. However, 
new features appear, which are due to very weak reso- 
nances of higher order in A. In particular, one clearly sees 
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FIG. 5: Growth rate j/A for the 2u) - resonance (O — 2cjofc) 
obtained with the stability analysis of Eq. I14II with so = 4, 
gn = 0.72E R and A < 0.01 (solid line) and A = 0.05 (dashed) 
and 0.1 (dot-dashed). 



narrow unstable regions at frequency loq and (u>o + u)i)/2. 
The growth rate and the width of these resonances are 
both of order A 2 , so that they are not accounted for by 
the two-mode approximation of section II I II The reso- 
nances at (u)q +u>i)/2 and 2loq exhibit a crossing, giving 
rise to a complex structure in the stability diagram. For 
our values of sq and gn, the crossing occurs at about 
k = 0.75qB- The effects of this crossing are also visible 
in the growth rate of the 2ujq resonance. In Fig. we 
plot j /A calculated at fi = 2<x>ofc for different values of 
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A. The solid line is the same as in the top-left panel of 
Fig.|21and corresponds to A < 0.01. The dashed and dot- 
dashed lines correspond to A = 0.05 and 0.1, respectively. 
The figure confirms that 7 remains directly proportional 
to A up to relatively large values of A, except in a nar- 
row region near k ~ 0.75<7s, in correspondence with the 
crossing point. 
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FIG. 6: Growth rate 7 as a function of the modulation fre- 
quency Q. The lines are the results of the stability analysis of 
sectionlTTlfor A — 0.2, so = 4 and gn = 0.72Er. In particular, 
the two solid lines are the resonances of order A at 2u>o and 
ljo + tt>i , while the dashed and dot-dashed lines are the weaker 
resonances of order A 2 at (luo+lji)/2 and (ujo+(j-'2)/2, respec- 
tively. Empty circles are extracted from the time evolution of 
the incoherent fraction SN/N in GP simulations on a cylindri- 
cal condensate, unbound along z and harmonically confined 
in the transverse direction. Solid circles come from GP sim- 
ulations as well, but for a condensate which is harmonically 
confined also along z, similar to the trapped condensates of 
the experiments of Ref. (lEl Htl Il7) . 



(cujk + ojj'k)T J ^_(k). Hence, the fact that for the 2uj q - 
resonance the growth rate 7 has a shape similar to that 
of utQk implies that T°^_ varies little with k. This is con- 
firmed by a direct calculation of T°^_ through Eq. 119fl . 
The result is shown as the solid line in the lower panel 
of Fig. [7| Conversely, the growth rates of the other reso- 
nances in Fig. [3] do not exhibit this behavior, indicating 
that the corresponding TV- have a nontrivial dependence 
on k. 

0.07 1 . 1 
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The growth rate can also be plotted as a function of 
the modulation frequency fi, as done in Fig. [5] where we 
show the results for A = 0.2. The solid lines represent the 
growth rate of the two resonances, of order A, at 2o>o and 
luq+lui, while the dashed and dot-dashed lines correspond 
to the weaker resonances, of order A 2 , at (uj + w\)/2 
and (cjq +£t»2)/2, respectively. The growth rate of the ujq 
resonance is very small on this scale and it is not shown. 
The points with error bars come from GP simulations 
and will be discussed in the next section. 

Coming back to the ^4^0 limit, we notice that, as 
shown in Fig. [21 the growth rates of the various types 
of resonances significantly differ in their dependence on 
the quasimomentum k. Apart from a slight upward bend 
at around half the Brillouin zone, the curve for the 2luq- 
resonance mirrors the form of the lowest Bogoliubov band 
(see Fig. nj. In the other cases, the fc-dependence of the 
growth rates is not that easily characterized. Some in- 
sight is gained from the two-mode approximation where, 
according to expression (|33|l . the fc-dependence of the 
growth rate on resonance is determined by the product 



FIG. 7: Coupling T + _ between counter-propagating modes of 
the first Bogoliubov band as a function of quasimomentum k 
at gn — O.IEr (top panel) and gn = 0.72Er (bottom panel) 
for different values of the lattice depth so = 4, 10 and 16. 
Lines: definition 1191 with the quasiparticle amplitudes taken 
from the numerical solution of the Bogoliubov equations JSJ . 
Points: Tight binding result 1421 . 

Up to now, our discussion has dealt with the proper- 
ties of parametric instabilities at a lattice depth which is 
modulated around so = 4. As so is increased, the width 
of the first Bogoliubov band decreases and, accordingly, 
the 2loq resonance is shifted to lower modulation frequen- 
cies f2. Moreover, an increase of the lattice depth changes 
the growth rates of parametric amplification. This can 
be seen in the two-mode approximation where the 2a>o _ 
growth rate is proportional to sq, the Bogoliubov fre- 
quency ujq and the coupling T~P_. In the lower panel of 
Fig. the quantity r°°_ is plotted for three different val- 
ues of s. As one can see, its s dependence is quite weak. 
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The product sqcoq remains also of the same order, since 
an increase of Sq is counterbalanced by a decrease of ujq . 
At zone boundary and for gn — 0, their product takes 
the values 4.2, 5.8, §.2Er for so — 4, 10, 16, respectively. 
This implies that the effects of a variation of so on the 
timescale of parametric instability of the 2u>q branch are 
not dramatic in this range of sq- For sufficiently large 
So one enters the tight binding regime, where loo is pro- 
portional to (see Eq.(|3HI)) and one gets an expo- 
nential decrease which eventually dominates the lattice 
depth dependence of the product so^o and hence of the 
growth rate. 

In section llVl above, we derived an expression for the 
coupling r°°_ in the tight binding regime (see Eas. (|4*2*|) - 
(I43|I % 1. In Fig. [7| we compare the tight binding results 
with those obtained from the numerical solution of the 
Bogoliubov equations © for different values of the lat- 
tice depth so and two different values of the interaction 
parameter gn. The numerical results (lines) approach 
the tight binding results (points) as the lattice is made 
deeper. For a given so the agreement is better when the 
interaction parameter gn is lower. For gn — O.IEr the 
agreement is reasonably good already at a lattice depth 
of s = 4, while deeper lattices are necessary to achieve 
the same degree of agreement for gn = 0.72.Er. This 
is consistent with general arguments on the applicabil- 
ity of the tight binding approximation [Isj ]. A repulsive 
mean-field interaction causes an effective smoothing of 
the periodic potential seen by the atoms and therefore 
the convergence to the tight binding regime is pushed to 
higher values so- 

With reference to Fig. 0we note that, at fixed so, the 
coupling T~P_ takes similar values for the two different 
interaction parameters gn. However, as pointed out in 
section III Bl the coupling between counter-propagating 
Bogoliubov modes, which is necessary for parametric am- 
plification, vanishes in the nonintcracting case. Numer- 
ically we find that the parameter gn has to be tuned to 
very small values (<C 10 -3 ) before a significant reduction 
of r5?_ is achieved. 



VI. COMPARISON WITH GP SIMULATIONS 

The above discussion of the parametric instability of 
the condensate is based on a suitable linearization of the 
ID GP equation @. We assumed that most of the atoms 
are described by an order parameter vE^o, which adiabat- 
ically follows the modulations of the lattice depth s(t), 
plus a small deviation S^. Under this hypothesis, we 
have found that <5M/ can exhibit an exponential growth 
as a consequence of the parametric instability of Bogoli- 
ubov excitations. The applicability of this approach to 
actual 3D condensates in experimentally feasible condi- 
tions can be tested by comparing its results with those 
of GP simulations, i.e., a direct numerical integration of 
the time-dependent GP equation (JIJ. 

For condensates subject to strong transverse confine- 



ment the integration of the GP equation can be efficiently 
performed by using the Non-Polynomial Schrodinger 
Equation (NPSE) introduced in Ref. 0. The key as- 
sumption is that the order parameter can be factorized, 
as in section [n] in the product of a Gaussian radial com- 
ponent and an axial wave function ^(z,i). Differently 
from section [HI using the NPSE one assumes the Gaus- 
sian to have a z-and i-dependent width, o~(z, f). The GP 
equation £[J thus yields 

2m \ a A a\ J 2-kg z 

(46) 

and a = ax(l + 2a7V|*| 2 ) 1 / 4 . In the geometry of the 
present work, differences between NPSE and the exact 
3D GP equation are negligible and solving the NPSE is 
much less time consuming. The effective ID GP equation 
(0) corresponds to the approximation a = a±. 

Let us first consider, as in section ^ a condensate 
which is unbound along z, subject to the external poten- 
tial 

V(z, t) = ^\r\ + s{t)E R sm 2 {q B z) , (47) 

with s(t) given in (|llfl . In order to compare the results 
with those of the effective ID GP equation J2J) we choose 
the parameters in (|46|l and (|47|l such that the condensate 
has the same interaction parameter gn, where n is the 
average linear density. 

We first calculate the ground state order parameter ^f 
in a static lattice of depth so and then we follow the evo- 
lution of \P under the modulation s(t), with given A and 
il. The Fourier transform of fy(z,t) gives the momen- 
tum distribution, n(k), as a function of time. The ini- 
tial momentum distribution is charaterized by the three 
peaks at k = and k = ±2qs, associated with the sta- 
tionary order parameter ^>q m the periodic lattice. As 
time goes on, we observe a relative oscillation of these 
peaks with the same frequency of the lattice modulation. 
This reflects the fact that \&o adiabatically follows the 
oscillations of the lattice depth. However, after a cer- 
tain time a parametric instability becomes visible, that 
is, the momentum distribution develops components at 
the momenta ±k of the amplified Bogoliubov modes. An 
example is shown in Fig. EI For even longer times, con- 
tributions from momenta throughout the first Brillouin 
zone start showing up due to higher-order harmonic gen- 
eration, Bragg reflection and nonlinear coupling between 
different modes (see 0] for a detailed discussion of the 
different stages of the evolution). 

The degree to which the parametric instability has 
evolved in time can be quantified by looking at the num- 
ber of particles SN contributing to the momentum dis- 
tribution away from the ^o components. In practice, 
we calculate the coherent part, N co h, by integrating the 
momentum distribution within the intervals ± A and 
±2qB ± A, with A = <7b/10, and then we define the 
incoherent fraction 8N/N = (N - N coh )/N. In Fig. 
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we plot 5N/N as a function of time as obtained with a 
modulation amplitude A = 0.2 for two different values 
of the modulation frequency f2 within the range of the 
2a>o-resonance. 
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FIG. 8: Momentum distribution of a condensate which is har- 
monically confined in the radial direction and periodic along 
z, obtained by performing GP simulations. The interaction 
parameter is gn = 0.72Er, where n is the average linear den- 
sity. The optical lattice has so = 4 and is modulated with 
A — 0.2 and = ujr. The modulation starts at t = and 
the momentum distribution is plotted at t — tj± (top panel) 
and t = ts (bottom panel), where A and B are the points 
indicated in Fig. |§] The arrows show the position of the para- 
metrically unstable modes. 

In order to observe the amplification, an initial seed is 
required. In the calculations a numerical noise is always 
present and is enough to trigger the instability with long 
modulation times. Adding an extra noise by hand in the 
initial order parameter anticipates the time at which the 
exponential growth becomes visible, but without affect- 
ing the growth rate. For the calculations in Figs. |S] and |§] 
the seed is added in the form of a white numerical noise 
in ty(z,0). This noise, which is not visible on the scale 
of Fig. HO is added in the interval — qs < k < qs and 
corresponds to 8N/N of the order of 10~ 4 . Only a small 
part of this noise is subsequently amplified, i.e., the one 
having the quasimomentum k of the resonant modes. By 
extrapolating the linear fit to t — we see that the seed 
in Fig. [21 is of the order of 10 -7 . 

We extract the growth rates 7 by plotting the data on 
a logarithmic scale and performing a linear fit in the time 
interval where the exponential amplification is visible, as 
indicated by the dashed lines in Fig. [5] The extracted val- 
ues of 7 can be directly compared to the results discussed 
in the previous sections, which were calculated by an ap- 
propriate linearization of the effective ID GP equation 
The comparison is shown in Fig. El The results of 
the numerical integration of Eq. Q are shown as empty 
squares. The error bars come from the uncertainties in 



the fitting procedure on the growth of SN/N in the sim- 
ulations. The overall qualitative agreement is good. The 
growth rate increases almost linearly with O, following 
the 2u>o resonance up to the band edge at k = qs- Then 
it suddenly drops. Weak A 2 -resonances become visible 
in the gap between the 2luq and the (loq + toi) branches. 
The agreement is especially good for the lowest 2uiq reso- 
nance. This confirms the accuracy, in this range of of 
the approximations made in section [H] namely, that the 
ID GP equation can be linearized in via the ansatz 
©, with a ^0 which adiabatically follows the modula- 
tion of the lattice, and that the inhomogeneous term in 
(|T2"j) is negligible. 
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FIG. 9: Incoherent fraction SN/N as a function of modulation 
time obtained by performing time- dependent GP simulations. 
Solid lines correspond to so = 4, A = 0.2 and two different 
modulation frequencies: f2 = lor and Q. — 2ujr. The growth 
rate 7 is extracted from a linear fit (dashed lines). The mo- 
mentum distribution of the condensate at the instants A and 
B is shown in Fig. [S] 

Now we consider a condensate which is also axially 
trapped. We solve again Eq. I|46fl. with the external po- 
tential 

V(z,t) = j(iolrl+Lu 2 z z 2 )+s(t)E R sm 2 (q B z). (48) 

We choose the parameters such as to simulate a sin- 
gle tube of Ref. 0. So, we use d — 413 nm, u> z = 
27rx84.6 Hz, uj± — 27rx36.5 KHz and the total number of 
particles N = 100. These values yield E R /h = 2ir x 3.34 
KHz. These are the same parameters used in our previ- 
ous work 01 • We extract the growth rate as done above 
for the simulation with the infinite cylinder. The results 
for the 2luq resonance are shown in Fig. El as solid circles. 
We notice that the number of particles per site in the 
trapped condensate is z-dependent, so that the compar- 
ison with the infinite cylinder of fixed linear density n is 
subject to some arbitrariness in the definition of the ax- 
ial average of n in the trap. Nevertheless the agreement 
is rather satisfactory and supports the physical pictures 
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discussed in 0] about the key role played by parametric 
resonances in the experimentally observed response to a 
lattice modulation in the superfluid regime 0, 0, • 
Finally, it is worth stressing that GP simulations ac- 
count for the effects of the nonlinear dynamics associated 
with the mean-field interaction term in the time depen- 
dent GP equation. In particular, when the population of 
a parametrically unstable mode becomes large, nonlinear 
processes yield a significant coupling with other nonres- 
onant modes and a consequent redistribution of energy 
over a wider range of quasimomenta, eventually causing 
a saturation of the incoherent fraction SN/N when it ap- 
proaches 1, as in the two curves plotted in Fig. [21 Effects 
of these processes are also visible in part B of Fig. 03 In 
the case of axially trapped condensates they are enhanced 
by the finite size of the systems, Z, which introduces a 
natural broadening of the order l/Z in the momentum 
distribution (see discussion in ^H). However, an impor- 
tant result of the present analysis is that such nonlinear 
processes do not prevent the parametric resonances to 
be observable, either by measuring the final effects of the 
amplification process (as, for instance, the total energy 
transferred to the condensate as a function of f2 and A as 
in 0,0,01) or by looking directly at the exponential 
growth of the unstable modes in the momentum distribu- 
tion. In the latter case, the identification of a peak with 
a rapidly growing height in the measured momentum dis- 
tribution would also allow one to use the parametric res- 
onances as a spectroscopy tool, namely to measure the 
dispersion relation ujjk in a parameter space where other 
spectroscopic techniques, such as two-photon Bragg scat- 
tering, are not easily applicable. The timescale and the 
size of SN/N predicted in this work seem well compatible 
with feasible experiments. 



VII. REMARKS ON THERMAL AND 
QUANTUM SEED 

As we already noticed, at the mean-field level the para- 
metric resonances occur only if the relevant Bogoliubov 
modes are initially occupied. Within the two-mode ap- 
proximation, discussed in section lTTTl we have shown that 
a special role is played by the effective seed rj defined in 
l|37p. This quantity is exponentially amplified in time. 

Different mechanisms can cause a non-vanishing rj. 
The condensate can be simply out of equilibrium at the 
beginning of the modulation due to non-adiabatic pro- 
cesses when it is loaded in the trap and in the optical 
lattices. This yields some excitations over the ground- 
state. A similar and controllable seed could be added 
to the condensate by exciting certain Bogoliubov modes 
by means of an appropriate external perturbation (e.g., 
two-photon Bragg scattering) before starting the lattice 
modulation. Both these types of seed can be described 
with the GP equation. On the contrary, the GP equation 
does not account for thermal and quantum fluctuations, 
which have to be treated with different methods beyond 



mean-field. 

In sections ITU and InD we treated the collective excita- 
tions of the condensate by linearizing the GP equation. 
We have also used the Bogoliubov quasiparticle ampli- 
tudes, solution of Eq. ©, as basis functions to represent 
the fluctuating component S^, the coefficients of the ex- 
pansion (|15|l being c-numbers related to the population of 
the collective modes. This allowed us to interpret the dy- 
namics in terms of coupling and growth of Bogoliubov ex- 
citations. A natural generalization consists in replacing 
these c-numbers with the quasiparticle creation and anni- 
hilation operators, as in the standard Bogoliubov theory 
[29l l3(| . In such a way one can treat both thermal and 
quantum fluctuations, under the only assumption that 
the total number of particles N is macroscopic and the 
fraction of noncondensed particles is much smaller than 
N. A direct connection between this quantum treatment 
of the operator S^ and our previous approach of section 
IIIII can be found by using the Wigner representation of 
the quantum fields 0, . In this way, the dynamics is 
still governed by the classical equations of sections ITU and 
IIIII but the condensate depletion is now included via a 
stochastic distribution of the coefficients Cjk- Exact re- 
sults can be obtained by averaging over many different 
realizations of the condensate in the same equilibrium 
conditions. In particular, one has 



fefc(0)) = , 
(c jfe (0)c*, fe ,(0)) = 



(c jk (0) Cj , k (0)) = 

1 1 

2 + e /31"jfc - 1 



(49) 

5jj'Skk> , (50) 



where (i = (fc s T) _1 . The 1/2 term in J5D|) describes 
quantum fluctuations and is the leading one at small T, 
while the second term accounts for thermal fluctuations. 
One can make use of (|49|) - (|5UII to evaluate the mean ef- 
fective seed defined in (|^7jl . finding 



(\v\ 



1 



1 



1 



(51) 



Equation (|51|l provides a way to extract information 
about the initial seed from the measurement of the expo- 
nential growth of the incoherent fraction. For example, 
within the two mode approximation described in section 
IIIII one finds that the number of particles in the modes 
j and j' at t ^> 7~a X , averaged over many realizations, 
grows as 



(SN) = C(\r, 



(52) 



withC = /^ 2 dz(|^ fe | 2 + |^| 2 + |^ fc | 2 + |^ fc | 2 ). So, if 
one is able to measure (5N) at different times during the 
parametric amplification and extract its value at t = 
from a fit, then Eqs. l|51[l and (jol'l) give direct information 
on the quantum or thermal nature of the seed, provided 
the Bogoliubov spectrum of the condensate is known. 

The thermal regime is obtained when kgT 3> hujjk 
and corresponds to the "classical" limit where the oc- 
cupation numbers of the relevant modes are large. In 
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this regime, the parametric response of the system to the 
lattice modulation can provide a thermometry. Let us 
suppose, in fact, that an experiment is set up to observe 
the exponential growth at = 2cjofc, as in the simulation 
of Fig. and that the dispersion law wofc is also known 
(either from calculations or by direct measurements of 
n(k)). Equation (JCTJl predicts (|?7| 2 ) ~ (l/2)k B T/fiuj k 
for ksT fuvok, so that the knowledge of the spectrum 
and the average seed vs. k provides an estimate of the 
temperature T. If the fit works well this means that the 
gas of excitations is in thermal equilibrium. This is an 
important result in view of the fact that the temperature 
we are speaking about can be much smaller than the 
temperatures which can be measured with current tech- 
niques, based on the detection of the thermal cloud. In 
the absence of a visible thermal cloud, parametric reso- 
nances appear to be an interesting method to make small 
thermal fluctuation detectable via a selective amplifica- 
tion. 

The other important limiting case is that encountered 
when kgT <C TujJjk- in this condition the initial popu- 
lation of the unstable modes is mainly due to quantum 
fluctuations and (\rj\ 2 ) ~ 1/4. The modulation of the 
lattice thus acts on the quantum vacuum, squeezing it 
and tranferring energy to the system. So, even at T = 0, 
one can parametrically excite the condensate by amplifi- 
cation of the vacuum fluctuations and this is a signature 
of its underlying quantum nature. This regime is anal- 
ogous to the parametric down conversion, a well-known 
phenomenon in nonlinear quantum optics |32| . A con- 
ceptually interesting point is that the amplification of 
quantum fluctuations acts as a source of entangled quasi- 
particles, analoguous to parametric sources of entangled 
photon pairs in quantum optics |33j . A particularly ap- 
pealing case is that of a modulation frequency Q such 
as to resonantly excite two counter-propagating quasi- 
particles in different bands, j and j'. In this case, the 
quasiparticle pair is represented by the entangled state 

\9) <x\j,k)\j',-k) + \j,-k)\j',k) . (53) 

A quite similar situation has been recently discussed for 
the generation of branch-entangled polariton pairs in mi- 
crocavities through spontaneous interbranch parametric 
scattering [34|. It is finally worth mentioning that the 
creation of quasiparticle pairs out of vacuum fluctua- 
tions can also be viewed as a manifestation of a dynamic 
Casimir effect: the environment in which quasiparticles 
live is periodically modulated in time and this modula- 
tion transforms virtual quasiparticles into real quasipar- 
ticles m. 



VIII. CONCLUSIONS AND OUTLOOK 

In this work we have explored the origin and the ef- 
fects of parametric resonances in elongated Bose-Einstein 
condensates subject to a ID optical lattice whose inten- 
sity is periodically modulated in time. We have used 
a suitable linearization of the time-dependent GP equa- 
tion to calculate the stability diagram and the growth 
rate of unstable modes. We have shown that the main 
mechanism of instability is a coupling between pairs of 
counter-propagating Bogoliubov excitations. This cou- 
pling is caused by the modulation of the background in 
which the excitations live. This picture emerges quite 
clearly when a two-mode approximation is considered. In 
this case, one can derive semi-analytic results for the cou- 
pling between Bogoliubov modes in different bands and 
discuss the convergence to the tight binding regime. The 
results have been compared with those of time-dependent 
GP simulations, for both condensates which are infinite 
and trapped in the axial direction. The overall agreement 
supports the interpretation of our previous GP simula- 
tions |14| . which were aimed at explaining the response 
to the lattice modulation observed in the experiments of 
Refs. Ellll [13 

in the superfluid regime. Our analysis 
also suggests the possibility to perform new experiments 
in order to observe the effects of the parametric reso- 
nances in a more controllable way as a tool for a novel 
type of spectroscopy and, possibly, for a characterization 
of quantum and/or classical seeds in actual condensates. 
Along this line, our final discussion about the role of the 
seed and the applicability of the Wigner representation is 
more intended to be a suggestion for future work rather 
than a quantitative analysis. This topic certainly de- 
serves further investigations. 

Our formalism can be naturally generalized to different 
types of periodic modulations and different geometries. 
An interesting example is the periodic translation (shak- 
ing) of the optical lattice as in the recent experiment of 
Ref. [3(| . Another simple case is the modulation of the 
transverse trapping frequency in elongated condensates 
with or without optical lattice. Work in this direction is 
in progress |37| . 
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